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ABSTRACT 

We report results concerning the distribution and properties of galactic extinction 
at high galactic latitudes derived from stellar statistics using the Sloan Digital Sky 
Survey (SDSS). We use the classical Wolf diagram method to identify regions with 
extinction, and derive the extinction and the extinction law of the dust using all five 
SDSS spectral bands. We estimate the distance to the extinguishing medium using 
simple assumptions about the stellar populations in the line of sight. 

We report the identification of three extinguishing clouds, each a few tens of 
pc wide, producing 0.2-0.4 mag of g'-band extinction, located 1-2 kpc away or 
0.5-1 kpc above the galactic plane. All clouds exhibit grey extinction, i.e., almost 
wavelength-independent in the limited spectral range of the SDSS. We discuss the 
implication of this finding on general astrophysical questions. 

Subject headings: ISM: dust, extinction, stars, galaxies: clusters: individual: Virgo 
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1. 



Introduction 



The Milky Way (MW) extinction was studied since the beginning of the 20 century to 
determine the wavelength dependence of the galactic extinction (the galactic extinction law), and 
to map the extinction in the different directions of the Galaxy. Schultz & Wiemer (1975) showed 
that for a wide spectral range the MW extinction law is A\ ~ A -1 , thus the extinction curve can 
be parameterized by a single variable Ry, the total-to-selective extinction ratio, defined as 



where E(B — V) = Ab — Ay is the colour excess or the reddening. Ry varies in different 
directions of the MW (Mathis, 1990), especially toward dense dust clouds. The mean value 
< R v >= 3.14 ± 0.10 (Schultz & Wiemer, 1975) and the typical galactic extinction law, that 
behaves as A\ ~ A^ 1 , are only MW averages. Wavelength-independent extinction A\ ~ Const, 
which is produced by large dust particles with sizes a > A, is called 'grey' extinction usually. 
For grey extinction the reddening converges to zero, and therefore Ry 3> 3.1. Knowledge of the 
interstellar extinction measurements and mapping is important for stellar studies and also for 
extra-galactic astronomy. Extinction maps are needed to correct the galactic and extra-galactic 
observations for the MW extinction. 

The first photometric method for measuring the extinction was proposed by Wolf (1923) and 
uses "Wolf diagrams". The method is based on the assumption that dust located between the 
object and the observer reduces the brightness of the object, thus the apparent magnitude m\ 
of the object increases by A\ (the extinction in magnitude units). In presence of extinction, the 
distance modulus is: 



where d is the distance between the object and the observer. 

It is possible to derive A\ by comparing the cumulative star count distribution in an 
extinguished region with that of a dust-free reference field. The two distributions are assumed to 
be identical, but the dust extinction shifts the distribution of the extinguished region, as shown in 




(2) 
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Figure 1. Since the stellar distribution in the MW is not uniform, the distribution slopes of the 
two regions may change with the apparent magnitude, but will remain parallel. The star count 
distribution changes with the galactic location, thus the two regions must be sufficiently close to 
reduce these differences. 
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Fig. 1. — Schematic Wolf diagram: star count distribution of an extinguished region and of a 
reference field covering the same area, as a function of magnitude. The distributions are assumed to 
be the same, but due to the presence of the obscuring material, the distribution of the extinguished 
region is shifted by A magnitudes. The magnitudes ml, m2 and the splitting point SP are defined 
in Section 2.2. 

The "Wolf diagrams" method was widely applied for extinction measurements between the 
1930s and 1960s (e.g., Bok 1956). As recently as 1996 the method used star counts performed 
by eye (Andreazza k, Vilas-Boas 1996). The availability of digital catalogs with positions, 
magnitudes, and sometimes colors, eased significantly this tedious and error-prone task {e.g., 
Cambresy 1999a, b; Nikolic et al. 2001; Krause et al. 2003; Froebrich et al. 2005). Cambresy 
(1999a,b) used USNO catalog star counts and the Wolf diagrams method to derive an all-sky 
extinction map and detailed extinction maps of giant molecular clouds. Froebrich et al. (2005) 
used 2MASS cumulative star counts to obtain detailed extinction maps of the Galactic plane. 

The implementation of the method has two main problems. The first is the proper choice of 
the unextinguished reference field. The second is the stellar density and the luminocity function 
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variations, which originate from the galactic structure itself. The latter problem was solved by 
Cambresy (1999a) using an adaptive area algorithm, which changed the analyzed region area to 
keep the star counts constant. 

The use of this method is still limited due to difficulties in the automatization of the algorithm, 
the large amounts of CPU and manpower requirements, and the lack of a deep optical or UV 
photometric stellar surveys covering both the regions of high and low galactic extinction. One 
advantage of the Wolf method is that it allows the estimation of the distance to the extinguishing 
cloud as the typical distance to the stars at the magnitude corresponding to the break in the 
cumulative star count distribution of the extinguished region. The other advantage of this method 
is that it deals directly with the extinction and not with the reddening, and is therefore sensitive 
not only to "normal" dust, but also to dust that produces no reddening. 

Another photometric method developed in early 1990's deals with star colours instead of 
magnitudes. Since for "normal" dust the extinction oc A -1 , extinction affects different colours 
differently. The extinction can be extracted from the colour-colour diagram of the extinguished 
region, since the extinguished stars are shifted in the colour-colour diagram by the colour excesses. 
Implementations of this method were by Lada et al. (1994), Lombardi & Alves (2001) and Sale et 
al. (2008). 

The advantages of the colour-colour method are simplicity, robustness and full automation, 
while the disadvantages are a strong dependence of the results on statistical and systematic errors 
in the photometric data, and the assumption of a wavelength-dependent extinction. 

Dust and gas are well mixed in the MW; one can therefore measure the gas column density 
associated with a dust cloud and estimate the extinction using a gas-to-dust ratio. The main 
disadvantage of the gas-based method is that the quantities are connected to dust properties 
only through average ratios, and thus may differ for particular locations. Goodman et al. (2009) 
compared three different methods of extinction estimation for the Perseus star-forming region: 
mapping using near-IR colour excesses, dust far-IR emission, and mapping the integrated intensity 
of the 13 CO emission. While all the methods produced morphologically similar maps, the gas-based 
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extinction values deviated greatly from the two others. 

The obvious advantage of the gas-based method is the data reduction simplicity. The 
convenience of this method led to the production of the first all-sky galactic extinction map. 
Burstein & Heiles (1978) solved the major problems of the reddening biases and varying 
gas-to-dust ratio by applying a complicated algorithm to the HI spectral data. They combined 
21-cm measurements with galaxy counts in the direction of 49 globular clusters, 84 RR Lyrae 
stars, and two early-type stars, and produced an all-sky galactic extinction map that had 10-20% 
accuracy and 13 square degrees spatial resolution for regions with |6| > 10° and z distances larger 
than 300 pc (Burstein & Heiles 1982, 1984). From the early 1980's to the end of the 1990's the 
Burstein & Heiles maps were the only tool for the galactic extinction correction of galactic and 
extra-galactic observations. 

The last method uses the far infra-red (FIR) dust emission that is dominated by thermal 
emission from dust grains. Schlegel et al. (1998, hereafter SFD) combined IRAS and COBE/DIRBE 
measurements to create a relatively high resolution (based on IRAS) and well-calibrated (based 
on COBE) all-sky galactic extinction map. The authors used 100 and 240 micron COBE/DIRBE 
measurements to create a map of the galactic dust temperature. The combination of this map 
with IRAS high-resolution 100 micron intensity map provided an all-sky galactic reddening map 
with 16% accuracy and 1/8 x 1/8 square degrees spatial resolution. 

The SFD maps still remain the most powerful and commonly used tool for the galactic 
extinction correction (e.g., the foreground galactic extinction estimates in the NASA/IPAC 
Extragalactic Database), but cannot be used to determine the distance to a dust cloud since they 
provide only a cumulative effect of FIR dust emission. It is clear that the SFD maps tend to 
overestimate the extinction. Arce & Goodman (1999) pointed out that the SFD maps overestimate 
the reddening by a factor of 1.3-1.5 in regions of smooth extinction with Ay > 0.5 mag. Yasuda 
et al. (2007) found that this deviation is connected with a degeneracy between the 100 and 240 
micron measurements, which does not allow a proper black-body spectrum fitting to the data. 
Yahata et al. (2007) found that the FIR emission from galaxies affects the SFD reddening maps 
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that tend to overestimate the extinction in regions with background galaxies. Here we shall 
demonstrate that SFD maps can in some cases also underestimate the extinction. 

Most methods listed above are not sensitive to grey extinction, since they deal with reddening 
and not with the extinction, or since they use average MW ratios. The exceptions are the Wolf 
diagrams and the FIR emission maps. The SFD maps can determine the presence of grey dust, 
but do not distinguish between "normal" and grey dust. This leaves the Wolf diagrams as the 
only method useful for detecting both "normal" and grey dust. 

Large 1 jiva dust grains, which produce grey extinction in the visible part of the EM spectrum, 
are only a tiny fraction of the MW grain size distribution (Mathis et al. 1977; Kim et al. 1994). 
Large quantities of grey dust grains could be found only in places where physical processes 
selectively support the production or survival of large-size particles. 

Grey extinction was detected by Sitko et al. (1984), Dunkin & Crawford (1998), Barge & 
Viton (2003), Patriarchi et al. (2003) and Skorzynski et al. (2003) in circumstellar disks of some 
MW stars. Patriarchi et al. (2003) pointed out the connection of grey extinction with some local 
features in the Carina association. The detections of Patriarchi et al. (2003) and Skorzynski 
et al. (2003) were questioned by Mafz Apellaniz (2005), who argued that grey extinction was 
determined mistakenly due to an incorrect treatment of the Lutz-Kelker bias. The presence of 
grey dust in a circumstellar disk was modeled and explained by Men'shchikov et al. (1999) as a 
result of an initial phase of planetary system formation, when large particles coagulate from small 
grains. Grey extinction was also discovered in high-redshift GRB host galaxies by Savaglio et al. 
(2003), Savaglio & Fall (2004), Stratta et al. (2004), Stratta et al. (2005), Chen et al. (2006), Li 
et al. (2008) and Perley et al. (2008). Savaglio et al. (2003) proposed that small dust grains are 
preferentially destroyed by GRBs and only the large grains survive. 

Greyish extinction (weakly wavelength-dependent) was found also in supernova remnants by 
Todini & Ferrara (2001), Bianchi & Schneider (2007), Nozawa et al. (2008). The origin of the 
greyish extinction in supernova remnants is probably the same as in the GRB host galaxies: small 
dust grains being selectively destroyed by the supernova outburst. 
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The main purpose of our study was to detect dust clouds at high galactic latitude and 
evaluate their properties. The SFD map predicts that in the studied region the extinction is 
A g i < m .4. Based on this, one would not expect to find high-extinction dust clouds in this region. 
We used the Wolf diagram method, measuring the extinction directly from the stellar statistics. 
Comparing our results to the SFD maps, we can confirm or reject the SFD predictions in this 
area. We identified three rather compact regions that show "grey" or greyish extinction. 

We also intended to detect possible significant galactic extinction in the direction of galaxies 
of the Hubble Space Telescope Key Project (hereafter KP) that aimed to determine Ho to ±10% 
(Freedman et al., 2001). If the extinction in the direction of the KP galaxies would not have been 
properly estimated, this could have induced a systematic error in the H measurement. The KP 
galaxies were chosen so that the SFD extinction prediction in their directions would be small 
(A g < m .l). Our results showed no significant extinction for three KP galaxies (NGC 4321, NGC 
4548, NGC 4535) located in the region of study (Figure 2). 



2. Data and reductions 

We used data from the SLOAN Digital Sky Survey (SDSS) DR6 with technical details given 
in Adelman- McCarthy et al. (2008). They include a list of objects recognized as stars by the SDSS 
pipeline with their equatorial coordinates, the apparent magnitudes in five SDSS bands - u'g'r'i'z' 
and the photometric errors in each band. The SDSS limiting magnitude (u' lim = 22.0, g' Um = 22.2, 
r iim = 22.2, i' Hm = 21.3 and z' lim = 20.5) is defined as the magnitude of the faintest point source 
with a probability of detection repeatability of at least 95%. Considering this, we excluded from 
our study all the stars that do not satisfy the SDSS magnitude limits in at least one of the bands. 
The region studied here is plotted in Figure 2 and covers ~ 650 square degrees of solid angle in 
the Virgo and Coma Berenices constellations. The approximate borders of the area are a = U. h 2 
to a = I3. h 5 and S = 0° to S = 26°. 
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Fig. 2. — The g' band SFD extinction map for the Virgo cluster region plotted in equatorial 
coordinates. The grey solid line shows the border of our database. White circles indicate the 
locations of four galaxies from the HST Key Project. 

2.1. Tiling and the "Bull's eye" algorithm 

We mentioned in the Introduction specific details of the Wolf diagrams method and some 
shortcomings of this method. We use the method here and solve the problems of the reference 
field and of galactic structure variations by using galactic coordinates and applying an algorithm 
called "Bull's eye". 

The coordinates of all the stars in the original database were transformed to the galactic 
coordinate system, and the entire region was divided into one square degree cells. We rejected 433 
cells at the border of the region, since parts may actually be beyond the region defined in celestial 
coordinates, and retained 2632 cells, each with ~1000 to ~5000 stars. Figure 3 shows the sky 
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Fig. 3. — The g' band SFD extinction map of the Virgo cluster region displayed in galactic 
coordinates. The grey solid line shows the rejected cells at the border of the database. The 
sample region is located above the grey line. 



The selection of a reference field for each cell, and the averaging of star count differences 
caused by the galactic structure, were obtained using the "Bull's eye" algorithm. We defined the 
cumulative star count distribution in cells within a small area centered on the extinguished region, 
and compared it with that of cells within a region surrounding the selected area. The reason for 
choosing two concentric or quasi-concentric regions for the analyzed and for the reference field is 
to eliminate differences in stellar density and luminosity distribution in different cells due to a 
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change of viewing direction within the Milky Way. Changes in the projected stellar density and 
the luminosity distribution with (I, b) are expected because the MW stars are not distributed 
uniformly. By selecting a surrounding concentric comparison area around the extinguished region 
we average over such effects. 

We assumed that the cells within the surrounding region are not extinguished. This 
implies that we assumed a small and relatively compact extinguishing cloud, well-separated from 
its surroundings. In the worst case, if the cells within the comparison region would also be 
extinguished, we would find only the differential extinction between the target and the comparison 
region. The distributions were normalized with respect to the areas of the two regions, so that 
only projected stellar densities per square degree were compared. 

Our stellar selection program allows choosing particular cells from the database and defines 
the surrounding comparison cells. It is therefore possible to choose regions with arbitrary shapes, 
not necessarily having circular symmetry. All the reference fields assumed to be unextinguished in 
this study were chosen to be within two degrees from the analyzed regions. 



2.2. Determining extinction and distance 

A difference between the cumulative star count distribution of the reference field and that 
of the extinguished region causes a horizontal shift of the latter by the extinction A (see Figure 
1). In practice, the distributions begin to diverge from a location at mi. The slopes of the two 
distributions return to the same value at m2, and from that point on the cumulative distributions 
remain parallel. The average magnitude between the two points mi and mi is called here the 
"splitting point" (SP). The estimation of the extinction A can be performed by averaging the 
distance on the magnitude axis between the two graphs, as done by Wolf (1923). We defined the 
extinction for each magnitude bin starting at mi and averaged the obtained extinction values. 
The error in extinction was estimated as the standard deviation of the extinction values for each 
magnitude bin from the average value. To compare the cumulative star count distributions in two 
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regions we performed a x 2 Goodness-of-Fit test. 

We determined the distance to the extinguishing medium as the typical distance to the 
stars located at the SP. This was determined from the apparent magnitude m\ and the absolute 
magnitude M\ of each star. To find the distance for each star we used only the SDSS g'r'i' data. 
We transformed the (g' — r') and (r' — i') colours into (B — V) and (R c — I c ) respectively, and 
determined the apparent magnitude my using the Fukugita et al. (1996) relations. We used 
the HR diagram data for main sequence stars from Cox (2000) and a cubic spline interpolation 
to derive My from the stellar colours. The (g 1 — r 1 ) and (r' — i') colours yielded two distance 
estimates for each star, marked as -D( 9 - r ) an d -D( r -i)- We stress that these distances are valid only 
for main sequence stars, which represent ;> 90% of stars (Mengel et al. 1979; Juric et al. 2008). 

For further distance measurements we used only one distance estimate for each star, provided 
that the relative difference between the two estimates is small: 



\e(D)\ = 



< 0.1 (3) 



D (9-r) + D (r-i) 

Stars with larger e(D) estimates have great distance differences due to unusual colour relations, 
thus they probably are not main-sequence stars or may be binaries. The number of stars with 
large e(D) estimates is consistent with the number of non-main sequence stars (Mengel et al., 
1979; Juric et al., 2008). The stars that do not fit the SDSS limiting conditions (~ 20% of the 
total sample), the stars with large relative difference between the distance estimates(~ 9 — 10%) 
and the stars too red for the Fukugita et al. (1996) relations (~ 7%) were rejected from further 
processing. The rejeted stars make up to ~ 40% of the SDSS stars in each sky tile. 

We estimated the extinction and the distance to the extinguishing medium in an alternative 
way with a simple model of a thin uniform dust screen with extinction located at a distance 
d. We modeled the star count distribution in the direction of the extinguished region by applying 
this thin screen, with parameters (d, A\), to the star count distribution of the reference field. We 
assumed that every star more distant than d is extinguished, while all stars nearer than d retain 
their apparent magnitudes unchanged. We used the stellar distances derived previously. We 
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estimated the best-fit distance and extinction for each dust cloud (d rn i n , A rnin ) by minimization of 
the fit in the (d, A) parameter space. We used the downhill simplex algorithm to minimize the 
X 2 Goodness-of-Fit value between the distribution of the extinguished region and the modeled 
distribution. The errors of the fit are determined with the "constant \ 2 boundaries" method 
described in Chapter 15.6 of Press et al. (1992). 

3. Analysis and results 

We selected three regions A, B and C from the surveyed area where the SFD predictions were 
large and checked their Wolf diagrams. These regions are shown in Figure 4 and are centered 
on (l,b): (248°,67°)a, (260°,65°) B and (297°,73°) c respectively. The areas of the regions are 
Qa ~ 35.9, Qb ~ 3.2 and VLc ~ 5.7 square degrees respectively. The SFD predictions for regions A 
and B are A g i ~ 0. m 20, and A g i ~ 0. m 25 for region C. Examples of selected Wolf diagrams for the 
regions A, B and C are shown in Figure 5. The errors of star count distributions are defined as: 

AiV(m) = VN/^j , (4) 

where s/N is the Poisson error of star counts and j is the number of cells, which represents the 
area of the analyzed region or the reference field. 

The three regions do not show significant extinction, despite the SFD prediction. We can set 
an upper limit to the extinction in these regions as A g t < 0. m l. A number of explanations can be 
proposed to understand the discrepancy: 

• There may be extinction in these areas, but it is much smaller than the SFD prediction, 
thus the SFD map overestimates the extinction. This was mentioned by Arce & Goodman 
(1999), Yasuda et al. (2007) and Yahata et al. (2007). 

• The extinction source in these areas may be very distant, perhaps even extragalactic, and 
the SDSS does not reach deep enough to allow its detection using a stellar statistics analysis. 
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Fig. 4. — Regions found to show no extinction where the SFD map predicted significant extinction. 
Grey rectangles represent the one-square-degree cells of the database. White solid rectangles denote 
the cells of the tested regions, white dashed rectangles denote cells of the reference fields. 

We estimate the limiting distance of our study to be of order 10 kpc, as the typical distance 
of the stars at the SDSS limiting magnitude in the region of study. 

• The determination of the extinction in this study is performed through a comparison of the 
selected region and its surroundings. If the reference field itself is extinguished, no extinction 
can be detected using the Wolf diagrams method. We can probably rule out this possibility, 
since the SFD prediction for the reference fields of all three regions is A g > ~ 0. m l indicating 
that, from a thermal IR point of view, there may only be a small amount of dust in these 
regions. 

We also checked the locations of three HST Key Project galaxies that are included in the 
survey region and were chosen to be in locations with small extinction: A g > < 0. m l. We did not 
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Fig. 5. — Examples of Wolf diagrams for the regions A, B and C. Solid lines are for the star count 
cumulative distribution in the analyzed regions, dashed lines are for the star count cumulative 
distribution in the reference fields. The title of each panel shows the value of the reduced \ 2 GoFT 
between the distributions. Panel (i) shows the g' band Wolf diagram for the region A, panel (ii) 
shows the i' band Wolf diagram for the region B, panel (hi) shows the v! band Wolf diagram for 
the region C. The behavior of all the SDSS bands is similar. 

find significant amounts of extinction at their locations. 

Finally, we scanned the entire surveyed region in g' band using the "Bull's eye" method 
described above. We performed the scan cell-by-cell, taking the central cell as analyzed region 
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Fig. 6. — Extinguished regions found in this study. White solid rectangles denote cells of the 
extinguished regions, white dashed rectangles mark cells of the reference fields. 

and the surrounding cells within two degrees from the central cell as a reference field. We used 
the x 2 GoFT on the normalized star count cumulative distributions in order to detect cells 
with significant deficiency of faint stars comparatively to the reference field. We identified three 
locations showing a deficiency of faint stars. We examined these regions using the other SDSS 
bands data and found that they showed a significant deficiency of faint stars in all SDSS bands. 
We combined the data for the extinguished cells in each of the regions to obtain more accurate 
Wolf diagrams. The three extinguished regions are shown in Figure 6 and their Wolf diagrams are 
shown in Figure 7. The errors of star count distributions are defined by Equation (4). 
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For each region we calculated the extinction in each SDSS band and compared 
it with that obtained from the SFD prediction using the extinction ratios for the 
SDSS bands derived by Fukugita et al. (2004) and Yasuda et al. (2007), namely 
A u > : A g , : A r , : A v : A z > = 5.2 : 3.8 : 2.8 : 2.1 : 1.5. We used the accuracy of the SFD 
map, 16% of measured reddening value (Schlegel et al., 1998), to calculate an error of the predicted 
extinction. Below we present the results for each of the regions together with a short discussion. 



3.1. First extinguished region 
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Table 1: Extinction calculated for the first extinguished region vs. that predicted by the SFD map 



Band 


Direct 


Model 


SFD 




m 16±0. m 05 


U - iy -0.07 


0. m 14±0. m 02 


Ag, 


m 24±0. m 05 


n mo/i+0.06 
u - *±-0.06 


0. m 10 ±0. m 02 


A r i 
A v 
A z , 


m 26±0. m 09 
0. m 28±0. m 07 
m 29±0. m 06 


0-281HI 

n mq7+0.06 
u - a '-0.04 

0. m 50l°,;°| 


0. m 07±0. m 01 
0. m 05 ±0. m 01 
0. m 04±0. m 01 



The region is centered on (l,b) = (292°, 72°) or (a, 5) = {12. h 6, 9.°5) and covers ~1.5 
square degrees. Figure 8 shows the detailed SDF map of the region and demonstrates that the 
extinguished patch does not coincide with an enhanced FIR emission region. The Wolf diagrams 
indicate a g' band extinction of about 0. m 24. The specific values, derived directly from the Wolf 
diagrams or via the thin screen model, are given in Table 1. We found that the distance to the 
extinguishing medium is 890lgQ pc using the direct method, or 1510lgQ P c using the thin screen 
model. We note that these distance estimates are very different and discuss this in Section 4. 
Using the different extinction values from Table 1 we plotted the wavelength dependence of the 
extinction in the SDSS bands, which is shown in Figure 9. 



3.2. Second extinguished region 
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Table 2: Extinction calculated for the second extinguished region vs. that predicted by the SFD 
map 



Band 


Direct 


Model 


SFD 


Ay 

Ay 
A r i 

Ay 

A z , 


0. m 22 ±0. m 02 
m 26±0. m 02 
0. m 26±0. m 03 
m 25 ±0. m 03 
m 29 ± 0. m 08 


n m A 7+0.28 
U - ^'-0.15 

n 778/10 + 0.06 

u - 4Z -0.08 

n +0.06 

u - 4i -0.09 

n mo-i +0.08 
u - 6L -0.05 

n mqr+0.08 
u - o,J -0.09 


0. m 21 ±0. m 03 
0. m 15 ±0. m 03 
0. m ll ±0. m 02 
0. m 09 ±0. m 01 
0. m 06 ±0. m 01 



The region is centered on (l,b) = (316°, 80°) or (a, S) = (I3. h 0, 17.°5) and covers ~2.3 square 
degrees. Figure 10 shows the detailed SDF map of the region from which it appears that, in this 
case, the extinguishing medium coincides with a higher FIR emission patch. The Wolf diagrams 
indicate a g' band extinction of about 0. m 26. The specific values derived directly from the Wolf 
diagrams or with the thin screen model are given in Table 2. The distance to the extinguishing 
medium is 1600^0 pc using the direct method, or 1920^3™ pc using the thin screen model. 
Using the different extinction values from Table 2, we plotted the wavelength dependence of the 
extinction in the SDSS bands, which is shown in Figure 11. 



3.3. Third extinguished region 
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Table 3: Extinction calculated for the third extinguished region vs. that predicted by the SFD 
map 



Band 


Direct 


Model 


SFD 


K' 

A g , 
A r i 
A v 
A z , 


0. m 26±0. m 04 
m 29±0. m 05 
0. m 26±0. m 07 
m 27±0. m 07 
m 28 ± 0. m 05 


n mnc+0.10 
u - zo -0.06 

0™40l°; ° 

n ni/|i +0.08 
u - 4i -0.08 


0. m 12 ±0. m 02 
0. m 09 ±0. m 01 
0. m 06 ±0. m 01 
0. m 05 ±0. m 01 
0. m 03 ±0. m 01 



The region is centered on (l,b) = (276.°5, 70.°5) or (a, 5) = (12. h 3,9.°5) and covers ~1.0 
square degrees. Figure 12 shows the detailed SDF map of the region. Here, as in the first 
extinguished region, the FIR emission is at a local minimum. The Wolf diagrams indicate a g' 
band extinction of about 0. m 29. The specific values derived directly from the Wolf diagrams or 
using the thin screen model are given in Table 3. The distance to the extinguishing medium is 
8OOI5Q pc using the direct method, or 1320^ pc using the thin screen model. Here also the 
distance estimates are very different, as found for the first extinguished region. This is relegated 
to the discussion in Section 4. Using the different extinction values from Table 3, we plotted the 
wavelength dependence of the extinction in the SDSS bands, which is shown in Figure 13. 

4. Discussion and conclusions 

We presented here a method to derive MW extinction using SDSS stellar photometry analyzed 
with Wolf diagrams. We showed that this method can not only estimate the amount of extinction 
in all five bands, but can also estimate a distance to the extinguishing medium, if this is assumed 
to be a compact and well-defined dust cloud. 

We applied the method to high galactic latitude areas predicted by SFD to have significant 
extinction, but we found none. We confirmed that in the direction of three HST KP galaxies there 
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is no significant extinction. 

We scanned the entire ~650 square degree area and found three regions with measurable 
extinction. We determined extinction and distance to the extinguishing medium by two alternative 
methods: direct measurement using the Wolf diagrams, and fitting a thin uniform dust screen 
model. The direct measurements produced significantly smaller distances to the extinguishing 
sources in all the regions when compared with those from the model fitting, while the extinction 
values estimated by the two methods are comparable. We suggest that the distances estimated by 
the model fitting should be adopted, since in the direct measurements we excluded ~ 40% of the 
stars because they were very red or had non-main-sequence colour relations or were fainter than 
the SDSS limiting magnitudes, and also because of the somewhat uncertain estimation of the SP. 

We also checked the effect of possible biases on our estimates. The distances to the 
extinguishing medium in the three regions may be underestimated by less than 3% due to the 
Malmquist bias. The effects of tiling, substructure and extremely red stars rejection were found to 
be negligible. The binarity bias could, in principal, produce an error in the distance estimation, 
but the extinction measurement, as performed here, would not be affected. We did not treat the 
binarity bias in our study, since it is very complicated and is controversial. 

We checked the effect of different median longitudes and latitudes of the analyzed regions and 
the reference fields on the extinction estimation. The differences between the median longitudes 
and latitudes of the analyzed regions and the reference fields were of order 0°.2. The effect of 
different centers was found to be insignificant, since its contribution to the extinction estimation 
is about ~0.01 mag for all the bands in all the regions. 

We also performed the following 'sanity' check: we split the reference field into two smaller 
and concentric reference fields and compared the distributions of the analyzed region with that of 
the first reference field and those of the two reference fields. In regions A, B and C the reduced 
X 2 value was less than 2 for the comparison between the distributions of the analyzed region and 
of the first reference field and for the comparison between the distributions of the two reference 
fields. In regions 1, 2 and 3 the reduced \ 2 value was less than 2 for the comparison between 
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the two reference fields, while the reduced \ 2 value for the comparison between the distributions 
of the analyzed region and of the first reference field was larger than 10. This shows that the 
analyzed region differs from its surroundings, while the first reference field is very similar to the 
second reference field. 

We performed the other sanity check by breaking the extinguished region into two smaller 
areas and comparing their star count distributions to the distribution of the unchanged reference 
field. The distributions of the smaller regions were similar and the extinctions derived from them 
were consistent with those of the entire analyzed regions. 

The three extinguishing clouds found here are a few tens of pc wide, and thus are not 
associated with circumstellar disks. The clouds produce 0.2-0.4 mag of (/-band extinction, and are 
located 1-2 kpc away or 0.5-1 kpc above the galactic plane. The extinction measured in all three 
regions does not follow the accepted MW wavelength dependence. The trend of the extinction law 
from the direct calculation and from the model fitting remains rather constant with wavelength, 
indicating grey extinction from large dust grains. This is the first clear evidence of grey extinction 
in the MW interstellar medium at high galactic latitudes. 

The results for the second extinguished region differ qualitatively from the results for the first 
and the third region and the distances estimated by two methods for it are similar. In the second 
extinguished region the SFD map shows stronger emission than in the surroundings, while in the 
first and the third region the emission is weaker than in the surroundings. The extinction law 
estimated by the model fitting and presented in Figure 11 shows a rising trend, although this is 
weaker than the mean MW extinction law trend. The extinction in this region is probably not 
grey, but weakly wavelength-dependent (greyish). 

The distances and extinction curves for the first and the third regions are very similar. The 
angular distance between these regions is 15° and their galactic latitude is b « 70° (see Figure 6). 
Adopting a mean distance of 1400 pc, we estimate the projected distance between the clouds to 
be ~ 130 pc, one order of magnitude smaller than the distance from the observer to these regions. 
We conclude that the dust clouds in these regions may possibly have the same origin and belong 
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to the same dust cloud complex. 

All three extinguished regions found in this study show wavelenght-independent or weakly 
wavelength-dependent extinction, are located more than 1 kpc above the galactic plane, and all 
are in the third quadrant. It is tempting to connect them with other strange entities in this part 
of the sky, such as possible cannibalized dwarf galaxies, but with the information collected so far 
this would be speculative. 

However, our finding of significant extinction at high galactic latitude, and the determination 
that this is grey extinction, should at the least inject a note of caution in studies of high-latitude 
objects and extragalactic sources. Milky Way dust can affect observational results in subtle and 
unexpected manners. 
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Fig. 7. — Examples of Wolf diagrams for the regions 1, 2 and 3. The notation is as in Figure 
5. Panel (i) shows the i' band Wolf diagram for the region 1, panel (iii) shows the r' band Wolf 
diagram for the region 2, panel (v) shows the i' band Wolf diagram for the region 3. Panels (ii), 
(iv) and (vi) are the enlarged plots of panels (i), (iii) and (v) respectively, ml, m2 and SP denote 
respectively the magnitudes ml, m2 and the splitting point defined in Section 2.2. 
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Fig. 8. — SFD Eb-v reddening map for the first region. White rectangles denote the one-square 
degree cells of the first region, grey dashed lines denote the equatorial coordinate grid. The values 
of the Eb~v for each "pixel" are according to the grey scale at the right. 
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Fig. 9. — "Extinction law" for the first extinguished region predicted by the SFD map, calculated 
directly and calculated from model. 
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Fig. 10. — SFD Eb-v reddening map for the second region. White rectangles denote the one- 
square degree cells of the second region, grey dashed lines denote the equatorial coordinate grid. 
Eb-v scale at the right. 
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Fig. 11. — "Extinction law" for the second extinguished region predicted by the SFD map, 
calculated directly and calculated from model. 
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Fig. 12. — SFD Eb-v reddening map for the third region. White rectangles denote the one-square 
degree cells of the first region, grey dashed lines denote the equatorial coordinate grid. Eb-v scale 
at the right. 
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Fig. 13. — "Extinction laws" for the third extinguished region predicted by the SFD map, calculated 
directly and calculated from model. 



